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C^^ ■ Abstract 

Cn I We study the fine tuning in the parameter space of the semi-constrained NMSSM, 

t^ ' where most soft Susy breaking parameters are universal at the GUT scale. We discuss 

the dependence of the fine tuning on the soft Susy breaking parameters M1/2 and ttiq, 
and on the Higgs masses in NMSSM specific scenarios involving large singlet-doublet 
Higgs mixing or dominant Higgs-to-Higgs decays. Whereas these latter scenarios allow 
a priori for considerably less fine tuning than the constrained MSSM, the early LHC 
results rule out a large part of the parameter space of the semi-constrained NMSSM 



X 

C^ ' corresponding to low values of the fine tuning. 



1 Introduction 

The first motivation for supersymmetric extensions of the Standard Model (SM) stems from 
the sohition of the naturalness or fine tuning problem in the Higgs sector of the SM [1-5] 
(besides the unification of the gauge couplings and the possibility to explain dark matter): 
In the SM with an ultraviolet cutoff A much larger than the electroweak scale Mz-, the 
bare Higgs mass squared rn^ must satisfy roughly m^ — A^ ~ M|. Hence rn^ must be of 
the order A^, but must be finetuned relative to A^ with a precision of the order M|/A^. 
This fine tuning is enormous for A of the order of a GUT scale, but only a tuning of 0{1) 
is considered as natural. Within supersymmetric (Susy) extensions of the SM with Susy 
breaking terms of the order Msusy, the necessary tuning between the parameters is of the 
order of M|/M| , and hence independent from an ultraviolet cutoff A. 

First results of searches for Susy by the ATLAS and CMS collaborations at the LHC, 
based on ~ 1 fb~^ of data taken at 7 TeV center-of-mass energy, imply lower bounds on 
Susy breaking gluino and up/down squark masses in the 1 TeV range [6, 7] (for the latest 
publications, see the ATLAS and CMS notes on the web pages [8,9]). These bounds reduce 
the phenomenologically viable range of parameters in Susy extensions of the SM. However, 
an obviously interesting question is the impact of these negative results on the necessary 
tuning between the parameters in the remaining parameter region. 

It is well known that the non-observation of a Higgs boson at LEP [10] implies already a 
"little fine tuning problem" in the Minimal SupersjTiimetric Standard Model (MSSM) (see 
e.g. [11-14]), where the field content in the Higgs sector is as small as possible, but large 
radiative corrections are required in order to lift the mass of the lightest neutral CP-even 
Higgs boson above the lower LEP bound. Large radiative corrections require relatively large 
Susy breaking top squark masses compared to the electroweak scale. Via loop diagrams, 
large top squark masses lead to relatively large soft Susy breaking Higgs mass terms, which 
require some tuning among the parameters of the MSSM such that the Higgs vacuum 
expectation values (vevs) are of 0{Mz)i well below the scale of the Higgs mass terms. The 
required tuning among the parameters of the MSSM is typically of the order of a few %. 

The "little fine tuning problem" of the MSSM, originating from LEP constraints, is 
alleviated in certain regions of the parameter space of the Next-to-Minimal Supersymmetric 
Standard Model (NMSSM). The NMSSM is the simplest Susy extension of the SM with a 
scale invariant superpotential, i.e. where the only dimensionful parameters are the soft Susy 
breaking terms. No supersymmetric Higgs mass term // is required as in the MSSM, since 
it is generated dynamically by the vacuum expectation value of a gauge singlet superfield S 
(see [15,16] for recent reviews). Together with the neutral components of the two SU(2) 
doublet Higgs fields H^ and Hd of the MSSM, one finds three neutral CP-even and two 
CP-odd Higgs states in this model. 

The additional coupling A of 5* to Hu and Hd can lead to a larger mass uih of the 
SM-like neutral Higgs boson if, and to mixings of the physical CP-even Higgs bosons in 
terms of the weak eigenstates S*, Hu and Hd, implying reduced couplings of the physical 
eigenstates to the Z boson. Both phenomena make it easier to satisfy the LEP bounds [17], 
and allow to alleviate the little finetuning problem [18]. (See [19] for an evaluation of the 
upper bound on rriH of about 140 GeV if A is required to remain perturbative below the 
GUT scale; if A is allowed to be larger, uih can be even larger [20].) 
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Moreover, H can decay preferably into lighter NMSSM-specific singlet-like Higgs bosons. 
In this case LEP bounds on ttih are lower, and again the required fine tuning can be 
considerably smaller than in the MSSM [18,21-24]. Consequently it becomes important to 
study the impact of the recent bounds from the LHC on the fine tuning within the NMSSM, 
which is the purpose of this paper. 

The most frequently used quantitative measure A for fine tuning is the maximum of 
the logarithmic derivative of Mz with respect to all fundamental parameters pf^^ (if the 
fundamental Lagrangian is given at the GUT scale) [11,25-33]: 



A = Max{A,^^^} 



A 



GUT 



d\xi{Mz) 



d\n.{p, 



,GUTN 



(1) 



(See [34, 35] for alternatives; sometimes M"^ instead of Mz is used in the argument of 
the logarithm, which leads to an obvious additional factor of 2. Subsequently we prefer 
to study linear relations between all masses and couplings.) A depends on the point in 
parameter space and is, roughly speaking, inversely proportional to the required fine tuning 
(as discussed above) between the parameters pf^^- Hence, for a given point in parameter 
space, A should be as small as possible, preferably of 0{\). Preferred regions in the 
parameter space spanned by pf^'^ are those where A is minimal (denoted by A^in)- In 
practice, the value of A depends on the choice of independent fundamental parameters 
defining the model, and on the implementation of phenomenological constraints as the 
dark matter relic density, the anomalous magnetic moment of the muon etc.. 

Including WMAP constraints on the dark matter relic density (but leaving aside the top 
Yukawa coupling ht in the list oipf^"^), A has been studied recently within the constrained 
MSSM (cMSSM, with universal soft Susy breaking terms at the GUT scale) in [30-33]. 
First investigations of the impact of the early LHC results on the Susy parameter space 
in the cMSSM have been performed in [36-39], in the cMSSM with non-universal sfermion 
masses in [40], and within the general MSSM in [38,39,41]. 

Compared to alternative procedures as likelihood scans and/or Bayesian techniques 
(see [42,43] for studies within constrained versions of the NMSSM, and [44] for a recent 
discussion) a disadvantage of A is that it does not allow to marginalise (i.e. to integrate 
over) parts of the parameters which, in turn, would allow to determine "most likely" values 
for given quantities as masses of specific particles, given all present experimental constraints. 
(Clearly, such predictions for "most likely" masses seem to be of limited use; for instance, 
they would have failed miserably if applied to the SM-like Higgs mass in the pre-LEP era.) 

One can also leave aside the issue of quantitative fine tuning, content oneself with the 
fact that the fine tuning in Susy is always much smaller than in the SM with a large cutoff 
A, and determine "most likely" values for parameters exclusively from best fits to data from 
electroweak precision experiments. The impact of recent LHC bounds on such best fits has 
been studied recently in [43,45-49]. 

On the other hand, the constraints from recent or future LHC results on the quantitative 
fine tuning measure A can contribute to the discussion on the impact of LHC results on 
Supersymmetry in general. Hence we will compare these constraints within the parameter 
space of the semi-constrained sNMSSM (where singlet-specific soft terms are allowed to be 
non-universal) to the cMSSM, obtained from the sNMSSM in the limit A, k — )■ 0. Therefore 



we study the dependence of A on the universal soft Susy breaking parameters M1/2 and rriQ 
(gaugino and scalar masses, respectively) and on the gluino and up/down squark masses. 
In particular we investigate the relevance of NMSSM specific scenarios in the Higgs sector 
for fine tuning, as scenarios with large singlet /doublet mixing and scenarios with dominant 
h -^ AiAi decays [16,21-24,43]. 

In the next Section we define the model, the procedure for the determination of A, and 
discuss some specific properties of A in the NMSSM. Our results and their discussion are 
given in Section 3, and we conclude with a summary in Section 4. 

2 Fine Tuning in the NMSSM and the MSSM 

The NMSSM differs from the MSSM by the presence of the gauge singlet superfield S. The 
Higgs mass term ^HuHa in the superpotential Wmssm of the MSSM is replaced by the 
coupling A of Hu and Hd to S and a self-coupling nS^, hence the superpotential W^mssm 
is scale invariant in this simplest Zs-invariant version of the NMSSM: 

Wnmssm = ^SHu ■ Hd + —S + htHu ■ QT^ + hbHd ■ Q B'^ + KHd ■ Lr'i^ , (2) 

where we have confined ourselves to the Yukawa couplings of H^ and Hd to the quarks and 
leptons Q, Tr, Bji, L and r^j of the third generation; a sum over the three generations is 
implicitly assumed. (In (2), for the first and the last time, the fields denote superfields.) 
Once S assumes a vev s, the first term in Wnmssm generates an effective /^-term 

/ieflf = As . (3) 

The soft Susy-breaking terms consist of mass terms for the gaugino, Higgs and sfermion 
fields 



-C 



h.c. 



MiBB+M2 J2 W^Wa + Ms J2 G"Ga 

a=l a=l 

Co = mljHX + mj,jHd\^ + ml\S\^ + ml\Q''\ + m'^lT^ 



R\ 



+ mj,\BJ,\+mi\L'\+mi\T'^\, (4) 

as well as trilinear interactions between the sfermion and the Higgs fields, including the 
singlet field 

-£trii = {htAt Q-H^T'^ + hAh Hd ■QB% + KA^ Hd-Lr'^ 

+ XAxHu-HdS + ^KA^S^'^+h.c.. (5) 

(Again, an effective MSSM-like i?-parameter B^s = A\ + ks is generated.) 

All parameters in the above Lagrangian depend on the energy scale via the corresponding 
renormalization group (RG) equations, which account for the dominant radiative corrections 
involving large logarithms. In the constrained NMSSM, one imposes unification of the soft 



Susy-breaking gaugino masses, sfermion and Higgs masses as well as trilinear couplings at 
the grand unification (GUT) scale Mqut^ 



Ml = M2 = M3 = Ml/2 , 

rriH^ = rriHa = mq = itit = ms 



rriL 



m-r 



mo, 



(6) 



Since the singlet superfield could play a special role (its couplings to a hidden sector, respon- 
sible for supersymmetry breaking, could differ from the MSSM-like fields), we will allow for 
non-universal singlet-specific soft terms ms and A^, at the grand unification scale. This is 
the so-called semi-constrained NMSSM, denoted by sNMSSM subsequently. Including the 
top quark Yukawa coupling due to its influence on the RG equations, the Lagrangian of the 
sNMSSM depends on eight parameters pf^'^ at the GUT scale: 



P: 



GUT 



Ml/2, mo, Ao, A, K, ms, A^ and h 



t ■ 



(7) 



The calculation of Ap^""" defined in (1) proceeds in two steps: First, we compute the 
variations 



^Susy 



d\n{Mz 



(8) 



ain(pf^^ 

with respect to the parameters p^ "^-^ at the Susy scale (the Susy scale is defined to be of the 
order of the soft Susy breaking terms). Subsequently these variations are contracted with 
the Jacobian 

(91n(pj"^^ 



J i' 



d\n{pj^^) 



(9) 



which takes care of the renormalization group running of the parameters between the Susy 
and the GUT scales. Then we obtain 



GUT 



E^ 



Susy 



J, 



«J 



(10) 



The parameters p^^^^ at the Susy scale are those appearing in the (effective) Higgs 
potential, whose minimization determines the vevs v^ of H^, v^ of H^ and s of S. These 
vevs determine, in turn, the quantities 



Mz = — -{vt + v^), tan/3 = — and /icg = As . 

2 Vd 
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Including the dominant top quark/squark induced radiative corrections, the three min- 
imisation equations Ei are given by [16] 



El : ml^ + /ie\ + >^'vl + '^{vl - vj) - ^/iefr(A, + ns) + '-^ In (M^Jml^) = 
E2 : m^ + ^2^ + X'vl + ^-l^(yj - vl) - ^^,^(A, + Ks) = , 



Vd' 



E:i : 



m, 



l + KA^s + 2K,'^s'^ + X'^(vf, 



2\kvm - A^A, 







(12) 



where Mgtop denotes an average value of the top squark masses. (It is not necessary to be 
more precise here, in contrast to the radiative corrections to the physical Higgs masses.) It 
is straightforward to express the vevs w„, Vd and s in terms of M|, tan/3 and fi^a with the 
help of these equations. 

Hence the relevant parameters Pi^^^^ at the Susy scale are given by (leaving aside the 
electroweak gauge couplings gi and g2, as well as Mgtop inside the logarithm) 

jr,f"''^ = niH^, niHa, ml, Ax, A^, A, k, and ht . (13) 

In order to compute the variations A^^**-*' (see (8)) with respect to these parameters, we 
use 

= ^^^ = E iS^P^' + §^^Mz + j^6 tan /3 + |^<5/iefr (14) 

i dPi ^ dMz 9tan/3 dfi^s 

for j = 1, 2, 3. Since all partial derivatives of the equations Ej can be computed explicitely, 
the three equations (14) can be solved for SMz (and, separately, for Sta.nf3 and 6fj,es) as 
function of all Sp^^^^ , which allows to determine the variations A.^^^^ in (8). 

At this stage it is useful to recall the origin of the "little fine tuning problem" in the 
MSSM. Neglecting the radiative corrections, the minimisation equations (12) of the Higgs 
potential imply, with fi^s = /i in the MSSM, 

In the absence of fine tuning, all terms on the right hand side of (15) should be of comparable 
magnitude, and no large cancellations should occur; hence both /i^ and |?Ti|^^| should not be 
much larger than (9(M|). However, from the RG equations one typically obtains fn^Hu ~ 
— Mg^^-Qp, which is often required to be much larger (in absolute value) than M|: At least 
within the MSSM, the SM-like Higgs scalar mass increases proportionally to In {M^^^^/m^^pj 
due to top/stop induced radiative corrections. Then, large values for Mstop are unavoidable 
in order to satisfy the LEP bound. Albeit large stop masses are consistent with the non- 
observation of top squarks, they would generate an uncomfortably large value for —Tn\j^ 
which has to be cancelled by yU.^ in (15). 

For large \Tn\j \ ~ /i^ one finds for tan^ /3 ^ 1, following (8) with i = niu^ or i = fi. 



I 2 I 9 

^susy ^ 2 ^" ~ A^"'y ~ 2-^ de) 



Accordingly large values for A^"*^-^ (leading, generally, to large values for Af^'^) reflect well 
the necessary fine tuning if |?Ti|^^| and hence /i^ are large. 

In the NMSSM /i is replaced by fics = As. For large |?ti|^^| ~ /i^^, the above reasoning 
remains essentially unchanged: For s ^ Mz (valid in most of the parameter space), E^ in 
(12) gives 

1 



^ ( -^« - \/Al -8ml]. (17) 



Replacing /i^ = /i^g and (17) for s in (15), one finds again from (8) 

I 2 I 2 

in the NMSSM. Hence, quite obviously, large values for |?7i|^^| are unnatural as well. How- 
ever, due to the NMSSM specific contributions to the Higgs masses and mixings or NMSSM 
specific Higgs decays, LEP bounds on the Higgs sector can be satisfied for smaller top/stop 
induced radiative corrections, hence for smaller values of Mgtop, allowing for smaller values 
for \mjjj and/i^g. 

It remains to express the variations of the parameters at the Susy scale in terms of 
variations of the parameters at the GUT scale, i.e. to compute the Jacobian Jij in (9), via 
the integration of the RG equations for the parameters. In the cMSSM (with boundary 
conditions at the GUT scale as in (6)) one can always write 

ml^ = a^^^ml + a^^^M^^^ + a^^^Al + a^^^Mi/a^o (19) 

and 

fx'' = bfxl, (20) 

where the coefficients a*^*) and b depend on the gauge and Yukawa couplings. 

In the typical case where all a^*^ in (19) satisfy \a^^^ I < 1? one can verify that all variations 
d\ia.{mH^)/d\n{pf^^), with pf^^ = mo, M1/2, Aq, are less than 1. At first sight, due to 
'-^Pi^'^ < ^m„^ for these parameters pf^^, this seems to reduce the necessary fine tuning 
in the MSSM. However, due to (91n(yU,^)/91n(/iQ) ~ 1, A^"*^^ ^ ^^^^ remains always large. 
Moreover, hf^'^ should generally be included in the list of parameters pf^^ [27,50], and the 
corresponding variation A^^""" can be large. This holds particularly in the so-called focus 
point region of the MSSM where mo ^ M1/2, Aq, and \a^^'^\ <^ 1 in (19) for specific values 
of hf^'^ (but a large derivative of o^^-'ttiq with respect to hf^^). 

In the sNMSSM, additional terms ~ ?7i| and ~ A^ appear on the right hand side of 
(19), which have little impact in practice. Instead of (20), the RG equations for A, k 
and m| will now play some role since, replacing (17) for s in fi^s, f^cs depends on these 
parameters (apart from a dependence on A^,). In fact one can verify that, in the MSSM 
limit of the NMSSM where A, k ^ 1, HeS satisfies the same RG equation as fj,. All in all 
we cannot expect dramatic effects on the fine tuning from the somewhat different running 
of the parameters between the Susy and the GUT scale in the NMSSM. 

In practice our computation of the different variations Ap^""" in the space of parameters 
P?^^ (7), in order to find its maximum A as function of i (see (1)) at a specific point in 
the parameter space, is performed as follows: For each such point, the code NMSPEC [51] 
inside NMSSMTOOLs [52,53] is used in order to compute the Higgs and sparticle spectrum 
including radiative corrections as described in these references. Constraints from LEP, B- 
physics and the anomalous magnetic moment of the muon are taken care of according to the 
latest updates given on the web site http://www.th.u-psud.fr/NMHDECAY/nmssmtools.html. 
(No constraints on the dark matter relic density are imposed; however, in most cases these 
could be satisfied by giving up the bino mass unification Mi = M1/2 at the GUT scale, i.e. 
chosing an appropriate mass for the lightest neutralino-like Susy particle (LSP) without 
impact on the results relevant here.) 
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For phenomenologically acceptable points, A^ "^^ is computed from the three minimiza- 
tion equations Ei in (12), following the procedure described above. The Jacobian Jij, i.e. 
the variations of the parameters at the Susy scale in terms of variations of the parameters 
at the GUT scale, is computed numerically from the two loop RG equations. This allows 
to obtain the necessary quantities Ap^""", whose maximum with respect to all parameters 

3 Results for the cMSSM and the sNMSSM 

To start with, we apply our procedure to the cMSSM, which allows for comparisons with 
the sNMSSM and the available literature. The relevant parameters pp^"'"(cMSSM) at the 
GUT scale are 

pfUT(cMSSM) = Ml/2, mo, Aq, /io, ^o, and h . (21) 

As usual, /i and B (and hence /iq and Bq) are determined by Mz and tan/3, but they still 
contribute to the definition of A. Next we scan over a grid of values of M1/2 and mo. 
For each set of these values, we scan over Aq and tan/3. Keeping only sets of parameters 
consistent with all phenomenological constraints, we look for values of Aq and tan /3 which 
minimize A as defined in (1) for fixed M1/2 and rriQ. The resulting minimal values of A 
can be represented in the plane M1/2 — niQ, or in the plane Mgiujno — ''^isquark (where mgquark 
refers to squarks of the first generation) in Figs. 1. 
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Figure 1: The minimal fine tuning A in the cMSSM. For each set of M1/2 and mo, A is 
minimized with respect to Aq and tan/3. Left panel: A in the plane M1/2 — mo. Right 
panel: A in the plane Mgiumo — 'Tisquark- Bounds within specific cMSSM scenarios from 
ATLAS [6] are indicated as black lines, and from CMS [7] as red lines (see text). 



In order to guide the eye, we have indicated lower bounds from ATLAS and CMS notes 



on analyses of jets and missing E^, based on an integrated luminosity L*"* ^ 1 fb~^: In 
the plane M1/2 — niQ, lower bounds from ATLAS [6] (interpreted within the cMSSM with 
tan/3 = 10, Aq = 0) are shown as a black line, and lower bounds from CMS [7] are shown 
as a red line. In the plane Mgiuino — '^squark, lower bounds from ATLAS [6] are shown as 
as a black line. (The latter bounds from ATLAS are obtained in a simplified model where 
squarks decay only into quarks + a neutralino with a branching ratio of 100 %.) 

In the white regions in Figs. 1, phenomenological constraints cannot be satisfied for 
any values of Aq and tan/3: Either a stau would be the LSP (left hand side of the left 
panel), or a charged slepton, chargino, neutralino or a CP-even Higgs boson is excluded 
by LEP2/Tevatron (lower part of the left panel), or squarks are excluded by CDF/DO (left 
hand side of the right panel) or are theoretically unaccessible (right hand side of the right 
panel) . 

Note that, for M1/2 < 350 GeV and uiq < 700 GeV in the left panel, A decreases hardly 
with decreasing Susy breaking parameters M1/2 and mo; the minimal value of A in the 
pre-LHC allowed region is about ~ 33. In this region of the parameter space we observe 
the "little fine tuning problem" of the MSSM due to the LEP bound on the SM-like Higgs 
mass, which hardly depends on M1/2 and ttiq- The impact of the LEP bound becomes clear 
once we minimize the fine tuning A for fixed lightest Higgs mass rriH (without imposing 
LEP constraints on the Higgs sector) as a function not only of Ao and tan/3, but also of 
Ml/2 and tjiq. The result is shown in Fig. 2. 
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Figure 2: A as function of the SM-like Higgs mass rriH in the cMSSM, without imposing 
LEP constraints on the Higgs sector. 



We see the strong increase of A with rriH for rriH > 108 GeV. Accordingly the LEP 
constraint rriH ^114 GeV implies A > 33 in agreement with Figs. 1 (for low values of M1/2 
and mo). In fact, itih is always just above 114 GeV in the entire planes in Fig. 1, once Aq 
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and tan/3 are chosen such that A is minimized for fixed M1/2 and rriQ, but LEP constraints 
are apphed. 

Hence, for low values of M1/2 and niQ (in the region M1/2 < 350 GeV and niQ < 700 GeV), 
A is determined by the LEP constraints on mn (implying the little fine tuning problem of 
the MSSM) and not much affected by the lower bounds on gluino and squark masses from 
early LHC searches: In the region above the ATLAS/CMS bounds, the minimal value of 
A increases only to ~ 40 — 50. 

We also see in Fig. 2 that A does not decrease systematically with rriH, once lower 
bounds on sparticle masses from LEP and the Tevatron are imposed: In order to minimize 
A, preferred values of nin would be in the range 100 — 110 GeV which is excluded in the 
MSSM, but not in the NMSSM (see below). 

For larger values of M1/2 or rriQ, the origin of the required fine tuning is different: Here 
it is simply the fact that the weak scale (determined essentially by — 2(/i^ + mjj )) is small 
compared to the Susy breaking scale, which requires some tuning between the parameters. 
Since mn^ at the Susy scale is closely related to the squark masses, A increases rapidly 
with m squark (for rUsquark ^ 1 TcV) as it is visible in the right panel in Figs. 1. 

The fine tuning in the cMSSM has recently been analysed in [30-33]. The procedure 
and precision in these papers is similar to ours, except that constraints on the dark matter 
relic density are applied in [30-33], but contributions to A from the top Yukawa coupling 
ht are left aside. From Fig. 7d in [31] we find, once LEP constraints are applied, a minimal 
value of A ~ 70 for not too large values of mo- Given that A in [31] is twice as large as 
our A defined in (1), this coincides well with the left panel of Figs. 1 for moderate values 
of Ml/2- However, for mo > 800 GeV, A decreases to ~ 10 in [31], whereas A increases 
with mo in the left panel of Figs. 1. In fact, for larger values of ttiq, A is dominated by 
contributions from ht whose absence in [30-33] explains the different results for A in this 
region. 

Next we turn to the sNMSSM. In various regions of the parameter space of the sNMSSM, 
unconventional properties of the Higgs sector allow to alleviate the LEP bounds, lowering 
the minimal possible values of A. We found it interesting to study these lower bounds on A 
separately for different scenarios in the NMSSM Higgs sector (see also [18]), since these will 
have very different implications for future Higgs searches at the LHC. Hence we distinguish 
subsequently the following two scenarios: 

(1) The lightest CP even Higgs boson Hi has a large singlet component (H/S mixing). 
This implies a reduced coupling to the Z boson, and allows for Hi masses well below 
114 GeV [10]. 

(2) A CP-even Higgs boson H decays dominantly into a pair of lighter CP-odd bosons, 
H — !■ AA, allowing again for Hi masses well below 114 GeV [10]. (We omit the index 1 of 
Ai in the following.) 

The search for the minimal fine tuning A in each of these scenarios is performed similar 
to the procedure in the cMSSM: Again we scan over a grid of values of M1/2 and mo. Now, 
for each set of these values, we scan over Aq, tan/3, A and A^ [k and m^ are determined 



by Mz and tan /3, but included in the definition of A) using a Monte Carlo Markov Chain 
(MCMC) technique. Keeping only sets of parameters consistent with all phenomenological 
constraints, we look for values of Aq, tan/3, A and A^ which minimize A as defined in (1) for 
fixed Ml/2 and ttiq, allowing us to represent the resulting minimal values of A in the plane 
Ml/2 — "^0, or in the plane Mgiuino — J^igquark- (In order to distinguish the scenarios above, 
we require essentially BR{Hi — )■ bb) > 0.7 for scenario (1), but BR{Hi — )■ AA) > 0.2 for 
scenario (2).) 

In the scenario (1) (H/S mixing), the corresponding results for A are shown in Figs. 3. 
Now the constraint on the left hand side in the left panel from the absence of a stau LSP 
has disappeared, since a singlino-like neutralino can be the LSP. From here onwards, the 
bounds from ATLAS and CMS are indicative only, since the signals for supersymmetry in 
the NMSSM can be different notably in the case of a singlino-like LSP [54]. 
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Figure 3: The minimal fine tuning A, defined as in Figs. 1 for the cMSSM, in the sNMSSM 
scenario (1). 



Compared to the fine tuning in the cMSSM in Figs. 1, we see that A can be considerably 
smaller for not too large values of the Susy breaking parameters M1/2 and rriQ [18]: A large 
singlet component of Hi (in the range 0.8 — 0.85) allows for lighter Hi masses compatible 
with LEP constraints, which reduces the required fine tuning, see Fig. 4 below. The mass of 
the second mostly SM-like CP-even Higgs boson H2 is always just above 114 GeV. This is 
now easier to satisfy than in the MSSM, since the doublet /singlet mixing shifts the mass of 
the mostly doublet-like Higgs boson upwards. The values of the NMSSM-specific coupling 
A do not have to be large to this end; its value is always < 0.01. 

Respecting just the pre-LHC phenomenological constraints, the fine tuning measure A 
can be as small as 14 for low values of M1/2 and mo in this scenario, but this is precisely 
the region which is constrained by the first unsuccessful searches for Susy at the LHC. 
(As stated above, these constraints depend on the decay properties of the u/d-squarks 
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and gluinos. Additional bino — )■ neutralino decay processes can reduce EJp^^^ signature in 
all sparticle decay cascades. Applying nevertheless the bounds for the cMSSM scenarios 
studied by the ATLAS and CMS collaborations to the sNMSSM, we find that the smallest 
admissible value of A becomes ~ 44 for M1/2 ~ 400 GeV, ttlq ~ 600 GeV, similar to the 
smallest admissible value of A in the cMSSM.) 

For larger values of M1/2 and ttiq, the fine tuning within this scenario (1) of the sNMSSM 
becomes similar to the one within the cMSSM: As explained above, the origin of the fine 
tuning is now the smallness of the weak scale with respect to Msusy and not the LEP 
constraints on the Higgs mass; hence the possibility to alleviate the LEP constraints within 
the sNMSSM is less relevant. 

In Fig. 4 we show A as function of the mass niH^ of the dominantly singlet-like Higgs 
state, minimizing A as a function of ttiq, iVfi/2, Aq, tan/3 taking into account LEP con- 
straints from Higgs and sparticle searches. (The irregular structures originate from the 
LEP constraints on Hi — )■ 66, which lead to irregular upper bounds on the coupling of Hi 
to the Z boson as function of mu^-) Again, A does not decrease systematically with de- 
creasing mn^ given the phenomenological constraints on sparticle masses, but the behaviour 
is somewhat different from the dependence of A on the SM-like Higgs mass in the MSSM. 
Here we find that A is minimal for mu^ below ~ 100 GeV which coincides with the Higgs 
mass range where a 2.3 a excess in H ^ bb is observed at LEP 2 [10]. (This excess could be 
explained here since, although Hi is dominantly a gauge singlet, it still has a nonvanishing 
- but reduced - coupling to the Z boson, and will decay dominantly into bb.) 




Figure 4: A as function of the mass niHi of the dominantly singlet-like Higgs state in the 
mixing scenario (1) described above, including constraints from LEP. 
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Next we turn to the sNMSSM scenario (2) involving light pseudoscalars A, into which the 
SM-like CP-even Higgs boson H can decay. Again the required fine tuning is reduced [21,22], 
since LEP constraints allow for smaller H masses than 114 GeV. These constraints depend 
on niA'- 

a) For niA > 10.5 GeV, A will dominantly decay into a pair of bb quarks. This signature 
has been studied by the OPAL and DELPHI groups at LEP [55, 56] implying rriH > 105 — 
110 GeV if H has SM-like couplings to the Z boson [10]. Still, this lower bound on niH 
allows for lower values of A. 

b) For rriA < 10.5 GeV, A decays dominantly into a pair of r leptons. The signature 
H — )■ AA — !■ 4r has recently been re-analysed by the ALEPH group [57] for rriA < 12 GeV, 
implying again lower limits on rriH- However, for rriA in the range 9 GeV < rriA < 10.5 GeV, 
A can and will mix strongly with the CP-odd rjb{nS) states [58] which imphes a considerable 
reduction of the BR{A — )■ t~^t~) [59]. Hence bounds on i7 — >■ AA — >■ 4r hardly constrain 
TTiH in this case; now we find that the dominant constraints result from the lower limits on 
rriH depending on remaining sub-dominant BR{H — )■ bb) [10]. 
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Figure 5: The minimal fine tuning A, defined as in Figs. 1 for the cMSSM, in the sNMSSM 
scenario (2). 



Our results for the minimal fine tuning A in the sNMSSM scenario (2) are shown in 
Figs. 5 in the same planes as before. Similar to the sNMSSM scenario (1), respecting just 
the pre-LHC constraints, the fine tuning measure A can be as small as 9 for low values of 
Ml/2 and rriQ. Applying naively the bounds for the cMSSM scenarios studied by the ATLAS 
and CMS collaborations to the sNMSSM, we find that the smallest admissible value of A 
becomes ~ 39 for M1/2 ~ 375 GeV, mo ~ 700 GeV, hence below the smalles admissible 
value in the cMSSM. 

In the region of low A, the BR[H — )■ AA) is larger than 80 %. Again, the values of 
the NMSSM-specific coupling A do not have to be large in order to favour H — )■ AA decays 
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once these are kinematically allowed; the value of A varies between 0.02 and 0.16 in Figs. 5. 
In Fig. 6 we show the minimal value of A as function of the mass rriH of the now SM-like 
state H, imposing LEP constraints. Again the preferred value of rriH is not always as small 
as possible; now the minimal fine tuning is obtained for rriH in the range 100 — 105 GeV, 
and increases again strongly for larger values of uih- 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
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Figure 6: A as function of rriH in the scenario (2), once LEP constraints are imposed. 



It is also interesting to study the dependence of the fine tuning on niA- At first sight, 
A depends hardly on rriA within the sNMSSM, where A^ differs from Aq at the GUT scale: 
A variation of ttia from to ~ 50 GeV (still allowing for H — )■ AA decays) corresponds to 
a variation of \Af^\ from ~ 5 GeV to ~ 40 GeV at the GUT scale, which has practically 
no effect on the fine tuning measure A. However, different values for ttia correspond to 
different lower LEP bounds on niH] lower bounds mn, in turn, affect A as in Fig. 6. The 
resulting minimal values of A as function of ttla are shown in Fig. 7. 

The structures in Fig. 7 can be explained as follows: For rriA ^ H GeV, LEP bounds 
on H2 — )■ HiHi — )■ 46 (here with H2 ~ H, and Hi ~ A) apply, see Figs. 3 in [10]. 
The corresponding lower bounds on mu are somewhat weaker for ttla > 30 GeV than 
for 11 GeV < m^ < 30 GeV, leading to a lower fine tuning for niA > 30 GeV where 
"iT^H ^ 106 GeV. (For m^ ~ 51 GeV, the lower bound on niH has a dip implying a dip 
in A. For tua ~ 57 GeV, A increases since rriH > 2mA > 114 GeV is required for 
kinematic reasons; here this scenario is obviously not preferred.) For rriA ^ H GeV, strong 
lower bounds on tuh from ALEPH [57] on the BR{H — ?> AA — )■ 4r) apply at first sight. 
However, due to a reduced branching ratio for A — )■ 2r form A — rjb mixing [59] , the window 
9 GeV < ruA ^ H GeV is hardly affected by these constraints, allowing for mn to assume 
values for which A is minimal according to Fig. 6. (For mA ~ 5 GeV, parameters have 
to be tuned in order to satisfy constraints from Bs — )■ fi~^fi~.) Hence we find two distinct 
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Figure 7: A as function of rriA in the scenario (2), once LEP constraints are imposed. 



regions for niA where H — )■ AA decays allow for a considerable reduction of the fine tuning. 

4 Conclusions 

We have studied the amount of fine tuning in the parameter space of the semi-constrained 
NMSSM, and compared it to the cMSSM. Representing the minimal fine tuning in the 
plane M1/2 — ttlq allowed us to study the impact of early LHC results. First we verified 
quantitatively, to what extend the NMSSM-specific scenarios in the Higgs sector allow to 
alleviate the LEP constraints with respect to the MSSM. We found indeed that a consider- 
able reduction of the fine tuning is possible in scenarios where the lightest CP-even Higgs 
state is dominantly singlet like, or in scenarios where H — )■ AA decays are possible, notably 
for rrtA ~ (10 ± 1) GeV and 30 GeV < m^ < 55 GeV. 

If one applies naively the early LHC constraints on the cMSSM in the M1/2 — niQ-plane 
to the sNMSSM, the impact on the minimal fine tuning is stronger in the sNMSSM since 
these constraints affect in particular the region of smaller Susy breaking terms where the 
fine tuning in the sNMSSM can be relatively small. Still, the necessary fine tuning in the 
sNMSSM in scenarios where H — )■ AA decays are possible is smaller than in the cMSSM, 
but only if the Susy breaking terms are not too large (hence if sparticles are not too heavy). 
Otherwise, if the Susy breaking terms are larger, the fine tuning does not originate from 
LEP constraints on the Higgs sector, but from the smallness of the weak scale with respect 
to the Susy scale. This hardly depends on details of the Higgs sector of the Susy model, as 
long as the fundamental parameters are defined at the GUT scale and influence each other 
through the renormalization group running between the GUT and the weak scale. Hence the 
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same question should be re-analysed in the NMSSM with gauge mediated supersymmetry 
breaking [60], where the Susy breaking parameters can originate from much lower scales, 
and where the Susy breaking Higgs mass terms can differ considerably from squark mass 
terms. 
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